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Abstract 

The aim of this paper is to estabhsh the framework of the enclosure method for 
some class of inverse problems whose governing equations are given by parabolic 
equations with discontinuous coefficients. 

The framework is given by considering a concrete inverse initial boundary value 
problem for a parabolic equation with discontinuous coefficients. The problem is to 
extract information about the location and shape of unknown inclusions embedded 
in a known isotropic heat conductive body from a set of the input heat flux across 
the boundary of the body and output temperature on the same boundary. In the 
framework the original inverse problem is reduced to an inverse problem whose 
governing equation have a large parameter. A list of requirements which enables 
one to apply the enclosure method to the reduced inverse problem is given. 

Two new results which can be considered as the applications of the framework 
are given. In the first result the background conductive body is assumed to be 
homogeneous and a family of explicit complex exponential solutions are employed. 
Second an application of the framework to inclusions in an isotropic inhomogeneous 
heat conductive body is given. The main problem is the construction of the special 
solution of the governing equation with a large parameter for the background inho- 
mogeneous body required by the framework. It is shown that, introducing another 
parameter which is called the virtual slowness and making it sufficiently large, one 
can construct the required solution which yields an extraction formula of the convex 
hull of unknown inclusions in a known isotropic inhomogeneous conductive body. 
AMS: 35R30 
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1 Introduction 

The aim of this paper is to establish the framework of the enclosure method [9] for possible 
application to some class of inverse problems whose governing equations are given by 
parabolic equations with discontinuous coefficients. 

The framework is given by considering a concrete inverse initial boundary value problm 
for a parabolic equation with discontinuous coefficients. The problem is to extract infor- 
mation about the location and shape of unknown inclusions embedded in a known isotropic 
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heat conductive body from a set of the input heat flux across the boundary of the body 
and output temperature on the same boundary. 

Let n be a bounded domain of R", n = 2,3 with a smooth boundary. We denote the 
unit outward normal vectors to dfl by the symbol u. Let T be an arbitrary fixed positive 
number. 

Given / = f{x,t), {x,t) G SilxjO, T[ let u — Uf{x,t) be the solution of the initial 
boundary value problem for the parabolic equation: 

dtU-V ■ 'jVu = inl^x ]0, T[, 

7Vw-i/ = /on9Qx]0, T[, (1.1) 
u{x, 0) = in fl, 

where 7 = 7(x) = satisfies 

(Gl) for each i, j — 1, ■ ■ ■ ,n ^ij{x) is real and satisfies Jij{x) — Jij{x) G L°°{fl); 

(G2) there exists a positive constant C such that 7(0;)^ • ^ > C|^p for all ^ G R" and a. 
e. X e Q. 

See [3] for the notion of the weak solution. This paper is concerned with the extraction 
of information about "discontinuity" of 7 from u and 'j'Vu-i' on dQx]0, T[ for some / and 
an arbitrary fixed T < 00. However, we do not consider completely general 7. Instead we 
assume that there exists an open set D with a smooth boundary such that D (Z ft and 
j{x) a.e. X E fl\D coincides with a smooth positive function jo{x) of x G Q and satisfies 
one of the following two conditions: 

(Al) there exists a positive constant C such that —(7(0;) — 7o(a;)/n)C ' C ^ C'ICP ^01 all 
^ G R" and a.e. x e D; 

(A2) there exists a positive constant C such that (7(2;) — ^o{x)In)^ ■ C ^ for all 

^ G R'' and a.e. x e D. 

Write h{x) = 7(0;) — 7o(x)/„ a.e. x e D. 
We consider 

Inverse Problem 1.1. Fix a T > 0. Assume that both D and h are unknown and that 
7o is known. Extract information about the location and shape of D from a set of the 
pair of temperature Uf{x,t) and heat fiux f{x,t) for {x,t) G 9f2 x ]0, T[. 

The D is a model of the union of unknown inclusions where the heat conductivity 
is anisotropic, different from that of the surrounding inhomogeneous isotropic conductive 
medium. The problem is a mathematical formulation of a typical inverse problem in 
thermal imaging. Note that in [4] a uniqueness theorem with infinitely many f for Inverse 
Problem 1.1 has been established provided h has the form bin with a smooth function b 
on D. Thus the point is to give a concrete procedure or formula which yield information 
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about the location and shape of D. Note that when n — 1, there are some results: the 
procedure in [1, 2] with infinitely many f and the formula of Theorem 2.1 in [10] with a 

single f. 

In [14] we considered Inverse Problem 1.1 in the case when 70 = 1 and n = 3 and 
gave four extraction formulae of some information including the convex hull of D. In 
Subsection 1.3 we will reconsider three results of those which employ infinitely many /, 
from the view point of the framework given here. 

Note that in this paper we do not consider the single input case. Formulae with a 
single / see Theorem 1.1 in [14] and [15] for a cavity which is the extremal case 7 = in 
D. Those can be considered as some kind of extension of the enclosure method for elliptic 
equations with a single Cauchy data started in [8] to the parabolic equations. 



1.1 A reduction to an inverse boundary value problem with a 
parameter-A general framework 

Define 

rT 

Wf{x,T) = e ^^Uf{x,t)dt, X E fl, T > 0. 

J 



The w — Wf satisfies 

(V • 7V - r)w = e-^^Uf{x, T) in Q, 



7V^w -ly^ r e"^7(^, t)dt for x e dQ. 
Jo 



(1.2) 



This motivates a formulation of the reduced problem given below. 
Given F{x,t) and g{x,T) with let w = w{x,t) be the solution of 

( V • 7V - r)w = e-^^F{x, r) in 
7Vtu ■ u = g{x,r) on dQ. 



;i.3) 



Inverse problem 1.2. Assume that F{x,t), D and h are all unknown and that 70 is 
known. Extract information about the location and shape of D from a set of the pair of 
w{x,t) and g{x,T) for x G dfl, r > 0. 

For this problem we propose the following general framework which reduces the prob- 
lem to construct a family of special solutions of an equation coming from the background 
body. 

Theorem 1.1. Assume that: there exist constants Ci and /ii such that, as r — > 00 

||F(-,r)|U.(^) = 0(e^^V^^); (1.4) 
we have a family {vr) indexed with r > tq > of solutions of the equation 

V • 7oV^; - rt; = m Q (1.5) 
satisfying the conditions, for some constants 112, A*3; A*4; C2 and C3 

||V^;,|U2p) = 0(e^^V'^^), (1.6) 
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Let g = g{x,T) be a function of x G dfl having the form 

g = ^(^r)-fo—\an, 
where ^ satisfies the conditions, for constants /i and /i' 



(1.7) 
(1.8) 



liminf T''|^(r)| > 

r — >-oo 



and 



|*(r)| = C>(t'^). 
Let w be the solution of (1-3). IfT satisfies 

T>Ci + Cs- 2C2, 

then 



lim — log 



/ _ dVr 

gvr - w-io—- 



dS 



(1.9) 
(1.10) 

(1.11) 
(1.12) 



Remark 1.1. It follows from (1.6) and (1.7) that 

lim -log ||Vv^||l2p) 



This is the meaning of C2 which is uniquely determined by ||Vvr||L2(D) with all r >> tq. 
Remark 1.2. It follows from (1.7) and (1.8) that C3 > C2. 

Remark 1.3. Prom the proof one obtains also the order of the convergence of the formula 
(1.12): 



(gVr - w7o^ ) dS 
an \ ov 



C2 + O 



logr 



This is important for a suitable choice of r in the case when the data is noisy. 

Here we present an application of Theorem 1.1 to the case when 70 = 1. Let c > 
and T >To — c~^. Let uj-^ e 5"""^, n>2 and satisfy u • uj-^ — 0. Set 



CT \ U! + i \ 1 -— U! 



C^T 



z satisfies 

We observe that 



Z ■ Z — T. 



(A - T)e^-^ = 0. 
Thus we take the family {vr)T>To Theorem 1.1 



(1.13) 

(1.14) 
(1.15) 



Vt-{x) = v{x; z) = e^'' 
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Note that: 

• it follows from (1.15) that the function e'"^*v{x; z) of (x, t) satisfies the backward heat 
equation 

{l\ + dt){e'^'v{x-z)) = Q. 

• the absolute value of e^^*f (x; z) coincides with e~'^(*~'^^''^) and this is a solution of the 
wave equation with the propagation speed 1/c. By this reason we call this c the virtual 
slowness. See also [10] for this interpretation. 

Recall the support functions of D and Q: 

hoii^) = supx • a;, hQ{uj) — supx ■ cu. 
x&D sen 

We have Iid^u) < hQ,{uj) for aU oj since D dVt. 

Applying Theorem 1.1 to (1.2), we obtain the following corollary. 

Corollary 1.1. Assume that 70 = 1. Let f be the function of {x,t) G dflx]0, T[ having 
a parameter r > defined by the equation 



du 



{x)(fi{t), 



(1.16) 



where a real-valued function ip e 1/^(0, T) satisfying the condition: there exists // e R 
such that 



lim inf 

T XXl 



[ e-^^ip(t)dt 
Jo 



> 0. 



Let Uf — Uf{x, t) be the weak solution of (1-1) for f — f{x, t). IfT satisfies 

T >2c{hn{uj)-hD{uj)), 

then 



lim — log 

r— >oo 2r 



/ f 

JdCl Jo 



-Tt 



—VT-{x)f{x, t; t) + Uf{x, t)-^{x) 1 dtdS 



;i.i7) 



;i.i8) 



:i.i9) 



Remark 1.4. There is no restriction on the position of center of coordinates relative to 
fl and D which affects on the sign of h£,{u). 

Since — /iq(— cj) = mix^^x • oj < inf^g/) x • oj < hc^Q), we have 

hn{oj) - hoiuj) < hn{uj) + hn{-uj). 

Thus T in Corollary 1.1 can be arbitrary small by choosing a small known c in such a way 
that, for example, 

T 

c < 

2{hn{oj) + hni-oj)) 

since from this one gets (1.18). Therefore if one wants to estimate D from the direction 
UJ by the formula (1.19) and the 'size' of fl at the cu direction, that is the quantity 
haif^) + h^{—oj), is too large compared with T (this is the most difficult case), the virtual 
slowness c in (1.13) should be chosen very small. This is one of the two roles of virtual 
slowness. In the next subsection we give another role of virtual slowness. 
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Let us explain how to deduce (1.19) from Theorem 1.1. Comparing (1.2) with (1.3), 
one knows that the g and F in Theorem 1.1 have the form 

fT Qy 

g{x,T)^ e-^^(f{t)dt-^{x), X edfl, T> (1.20) 

and 

F{x,t) =Uf{x,T). (1.21) 

Thus (1.9) is satisfied with the same as (1.17); (1.18) is satisfied with /i' = —1. We 
have to know also an estimation of Ci in (1.4) from above. Prom [3] it follows that 

im-,r)|U2(n) = O(||/|U.(0,T;if-i/2(af2))) (1-22) 

and thus one can choose 

Ci = cha{uj). 

We have also 

C2 = choioj) < chn(uj) = C3. 
Since C1 + C3- 2C2 = 2c{hn{uj) - hniuj)), (1.11) becomes (1.18). 

Remark 1.5. Note that the choice (1.13) of z satisfying (1.14) goes back to [11] in which 
an application of the enclosure method for inverse source problems for the heat equations 
are given. The point of the choice is: the growing orders of and ^ • ^ as r — > 00 are 
same. See also [12] for an application to the so-called inverse heat conduction problem. 
In one-space dimensional case, in [10] instead of (1-13) the z having the form 

has been used. In this case Re (z-z) — r. Formula (1-19) can be considered as an extension 
of a result in one-space dimensional case [10]. It means that if one could always control 
the initial temperature in the process of all possible measurements at the boundary to be 
zero, then one can extract the convex hull of unknown inclusions. 

Since / is complex-valued, Uf on dflx ]0, T[ can not be directly measured and should 
be computed from real data via the formula 

^/ = ^Re/ + ^^Im /■ (1-23) 

This is a consequence of the zero initial data. Thus the zero initial data is essential 
for this procedure. Note also that since both Re / and Im / are highly oscillatory as 
r — > 00 with respect to the space variables, it will be difficult to prescribe those fiuxes 
on the boundary directly. Instead one has to make use of the principle of superposition to 
compute the right-hand side of (1-23) on dfl from experimental data which are generated 
by finite numbers of independent simpler input fluxes on dD,. Needles to say, for this 
procedure the zero initial data are also essential. 
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1.2 The case when 70 is not necessary constant. 

It is possible to to extend Corollary 1.1 to the case when 70 is not necessary a constant. 
For simphcity of description assume that 70 — 1 G C^(R"). We construct a special 
solution of the equation of (1.5) which has the form 

v{x) ^ 



as r — > 00, where z is given by (1.13). 

Following [19], we make use of the change of the dependent variable formula (the 
Liouville transform): 

^ V • 7oV ( ^ ■]=A-V, (1.24) 




/7o 

where 



We find the special solution of (1.5) having the form 

pX-Z 

[l + ez): 



V^—^. (1.25) 

VTo 




7o 

where is a new unknown function. It follows from (1.24) and (1-14) that the equation 
for e — becomes 

(1.26) 

Thus the problem is to construct a solution of equation (1.26) such that ^ in Q 
as T — > 00. In general this is not an easy task because of the growing factor r on the 
zeroth-order term in (1.26). However, we found that: if the virtual slowness c in (1.13) is 
sufficiently large and fixed, then one can construct such a solution for all large r >> 1 with 
an arbitrary small by using a combination of the Fourier transform and perturbation 
methods in [19] for the construction of the so-called complex geometrical optics solutions 
of the equation V ■ 70 Vf = 0. Its precise description is the following second result. 

Theorem 1.2. Let —1<6<0 and a,b E C^(R"). Given 1] > there exist positive 
constants Cj = Cj{a,b,fl,6,ri), j = 1,2 such that: if c > Ci and r > C2, then c^t > 1 
and there exists a unique G L^(R") with z given by (1.13) such that 

{A + 2Z-V -Ta-b)e;,^Ta + b inlC. (1.27) 

Moreover, e^jn can be identified with a function in C^{Q) and 

\\^z\\l'=°{q,) + ||Ve2||Loo(f2) < r]. (1-28) 

This theorem indicates the important role of the virtual slowness c when 70 is not 
necessary constant. It is not an accessary! To the best knowledge of the author this idea, 
that is, choosing a large c and fix, never been pointed out. 
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Having this theorem, we obtain a result which corresponds to Corollary 1.1. 
Let < 77 << 1 and fix a c > Ci in Theorem 1.2. Let a — (I/70 — 1) and b is given 
by (1.25). Let be the solution of (1.27) constructed in Theorem 1.2. Define 

Vr{x) = f (1 + e,{x)), xeQ, r > C2. (1.29) 



The function e '^*Vr{x) satisfies the backward parabolic equation 

(V-7oV + 9t)(e-^V(a;)) = 
and its absolute value has the form 

p—T{t—CX-U)) 

\l + eJx) 



This again supports the name virtual slowness of c. 

It is easy to see that the family (^Vr)T>C2 satisfies (1.5), (1.6) and (1.7) with C2 ~ 
choico) , (1.8) with C3 = chQ{u) and (1.4) with Ci = chn{uj). Thus applying Theorem 
1.1 to this case, we obtain the following corollary. 

Corollary 1.2. Assume that 70 — 1 G C^(R"). Fix the virtual slowness as c — Ci, where 
Ci is just the same as Theorem 1.2. Let f be the function of {x,t) e 9Jlx]0, T[ having a 
parameter r > defined by the equation 

dv 

f{x,t) = —{xMt), 

where v — Vr is given by (1.29) and a real-valued function </? e i^^(0, T) satisfying the 
condition (1.17) for a /j, & R. Let Uf — Uf{x,t) be the weak solution of (1-1) for f — 
f{x,t;T). If T satisfies 

T > 2c{ha{uj) - hoiuj)), (1.30) 



then 

1 r rT / /te- \ 



lim — log 



Lnlo ^ ^* ^-^^T(a;)/(a;,t;r) +M/(a;,t)7o^(a;)j (itd^ 



In Corollary 1.2 c = Ci and thus (1.30) should be considered as a restriction on the 
length of the time for data collection. A sufficient condition to ensure (1.30) is 

T >2Ci(hn{uj) + hn(-uj)). 

Note also that the center of coordinates in Theorem 1.2 and Corollary 1.2 is free from D 
and ft. 
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1.3 Real VS Complex 

Replace the conditions (1.4), (1.6), (1.7) and (1.8) with the following ones, respectively: 

||F(.,r)|U.(^)=0(e^^^r''^); (1.31) 

||Vt;,|U2(D) = 0(e^^^r^^); (1.32) 

W'^v^LHD) > C'e^'^T^^- (1.33) 

||^r||ffi(n) = C>(e^^^^T'^^); (1.34) 

Then instead of (1.12) we have, for any fixed T > without (1.11) and exactly same g 



satisfying (1.9) and (1.10) for some e R 

1 



lim ^ 

r— J-oo 2^yT 



log 



an 



9V7 



^70^) dS 



:i.35) 



Since the proof is simpler than that of Theorem 1.1, we omit its description as a theorem. 
Exphcit examples of v satisfying (1.32), (1.33) and (1.34) in the case when 70 = 1 and 
n — 2) are the following: 



v{x]T,p) 



e^^-^, a; e R^ we S^, 

p-\/t\x-p\ _ 

-,xe R^\M, peR^\n, 



\x — p\ 



v{x;T,y) = 



\x - y\ 



, X e R' \ {y}, v{y; r, y) = 2t, y e R^ 



These are all real- valued functions and not oscillatory as r — > 00. We think that this 
non oscillatory character is an advantage in computing the left-hand side of (1.12). 

If g and F are coming from (1.20) and (1.21) for / given by (1.16) for a cp satisfying 
(1.17) for a // e R, then one can choose 



Co 



hoioj), \i V — v{x;t,uj), 
-doip), iiv ^v(x;r,p), 
, Roiy), iiv = v{x;T,y), 



where d^:{p) and R^{y) denote the distance of p to * and minimum radius of the open ball 
that contains * and centered at y. See Lemma 3.1 in [14] and Proposition 3.2 in [13] for 
these facts. 

By virtue of (1.22) one can choose 



hn{uj), if V = v{x;t,uj), 
-da{p), ifv ^ v{x;t,p), 
. Rn(y), if = v(x;r,y). 
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Then, formula (1.35) reproduces Theorems 1.2-1.4 in [14]. 

However, when 70 is not necessary constant, to construct a suitable v one has to solve 
the eikonal equation 

7o Vf ■ Vv = 1 inQ 

and corresponding transport equations. However, this is not a simple matter in general 
because of the complicated behaviour of the characteristic curve x = x{t) under suitable 
initial conditions on {x{t),^{t)): 

-r = V 
dt 



Thus Corollary 1.2 suggests that when the background body is isotropic, however, not 
necessary homogeneous, the use of complex geometrical optics solutions are better than 
that of geometrical optics solutions. 



2 Proof of Theorem 1.1. 

In what follows for simplicity we wright Vr = v and 7o/„ = 70- 

Let R-yir) and Rioir) denote the Neumann-to-Dirichlet maps on dVt for the operators 
V • 7V — r and V • 70V — r, respectively. We have 

( dv \ _ 
R-fo{T) (7o^|an 1 ^v\d^, R^{T)g ^ p\dn, 

where p solves 

(V • 7V - t)p = infi, 

(2.1) 

7Vp • u — g{x, t) on 

Our starting point is the following identity which is an easy consequence of equations 
(1.3), (1.5) and integration by parts: 



+ / (7 ~ 7o) Vtj • V{w — p)dx + e""^^ / F{x, T)v{x)dx. 
Jn Jn 

Define e = w — p. It follows from (1.3) and (2.1) that e solves 

(V • 7V - T)e = e-^^F{x, r) in Q 

7Ve ■ 1/ — ondQ. 

Since r > 0, it is easy to see that 

||Ve|U2(f,) < Ce-^^T-V2||F(-,T)|U2(^) 



(2.2) 
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-r(^-Cl)^Ml-l/2^ 



and from (1.4) one gets 

\\^4LHn) = 0(e" --r" 
From this together with (1.4), (1.6) and (1.8) we obtain 

/ (7 - 7o) • Vedx + e"""^ / F(x, T)v{x)dx 
Jn Jn 

= 0(e"'^^'^"^i"^2V''i+''2"^/^) + 0(e"^('^"'^i "'^3 V''! +''*). 
A combination of (2.2) and (2.3) gives 

r / _ dv \ f dv ( dv \ 



(2.3) 



(2.4) 



The following type of estimates now are well known and it is a consequence of Propo- 
sition 2.1 in [14] which goes back to [7] and one of assumptions (Al) and (A2). See also 
[16] when 70 = 1 and h = {k — l)/„ with a positive constant k. 

Lemma 2.1. There exist C > and C > such that for all v satisfying (1-5) 



< 



f dv ( dv \ 



< C"||V^;||? 



(2.5) 



Note also that: 

• (1.9) is equivalent to the statement: there exists a Co > such that, as r — )■ 00 

Co<r^|^(r)|. (2.6) 
From the right-hand side of (2.5), (1.6), (2.4), (1.10) and (2.3) we have 

= 0(e^'^2r^2M2+M'(]^ _^ g-T(r-Ci+C2)^Mi-M2-;^'-l/2 _^ g-r(r-Ci+2C2-C3)^/ii+M4-2/i2-M'))_ 



One the other hand, from the left-hand side of (2.5), (2.6) and (1.7) we obtain 



(2.7) 



/ ( 



gv - w'yo^j dS 



_ (Jin^2C2T^2lJ.3-^^-^ _^ Q^g-T(T-Cl+C2)^W+/"2-2/i3+/i-l^ _^ Q^g-r(T-Cl+2C2-C3)^/il+M4-2/i3+/i^^^ 

(2.8) ' 
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where C" = CCo{C"f > 0. 

Now formula (1.12) is a consequence of (1-11), (2.7) and (2.8) provided 

T > max(Ci - C2, Ci + C3 - 2C2). (2.9) 

However, by Remark 1.2 we have (Ci + C3 — 2C2) — (Ci — C2) = C3 — C2 > and thus 
(2.9) is nothing but (1.11). 

□ 



3 Proof of Theorem 1.2 

3.1 A special fundamental solution 

Given F we construct a solution of the inhomogeneous modified Helmholtz equation 

(-A + T)^; + F = OinR" (3.1) 

that has the form 

v{x) = e"-"^(a;) (3.2) 

where the complex vector z is given by (1.13). 
Write 

F{x)^e^-'f{x). (3.3) 

Then if ^ satisfies the equation 

- - 2;^ • + / = in R", (3.4) 

then the v given by (3.2) satisfies (3.1) with F given by (3.3). 

We construct a solution of (3.4) by using a special fundamental solution of (3.4). Set 

a(0 = ier-2i^-e, eeR"- 

We have 

Reg,(0 

and 



- c't' ( 1 - 



1 



lmg,(0 = -2c™ 

Thus the set Q~^{0) = G R"^ \Qz{0 — 0} consists of the circle on the plane u ■ ^ = 



centered at — cry 1 — (c^r)~^C(; with radius cry 1 — (c^r)"^. Moreover VReQzi^) and 
VlmQz{^) on Q^^{0) are linearly independent. Thus —l/Qz{0 is locally integrable on 
R" and can be identified with a unique tempered distribution and its inverse Fourier 
transform 

is well defined. This g = Qz satisfies, in the sense of tempered distribution 

-A^ -2z-Vg + 5{x) = inR'* 
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and thus, in the sense of Schwartz distribution 

(- A + r) (e'^-'g,) + 6{x) ^ in R". 

The Gz{x) — e^'^gz{x) should be called the Faddeev-Green function. 

Prom [19], one knows that: by virtue of Theorem 7.1.27 of [5] given / G L5_|_i(R"') 
with — < 1 < 5 < the solution of equation (3.4) is unique in L^(R'^); if / is a rapidly 
decreasing function, then the unique solution of (3.4) is given by 

* = (3.6) 

its first and second derivatives belong to L^(R") and satisfy 

IP"*IIl|(R") < C«(^,5)||/|U2^^(Rn), |a| < 2. (3.7) 

Thus the operator :/ i — > ^ has the unique continuous extension as a bounded linear 
operator of L^_,_^(R") to L^(R"). We still denote the operator by the same symbol Qz* f 
which yields the unique solution of equation (3.4) in L^(R") for all f e Lj_^_-^(R^) . 

However, for our purpose, one has to clarify the behaviour of Ca{z,S) in (3.7) as 
r — )■ oo. This is not the well known case z ■ z = —k'^ with a A; > which appeared 
in inverse scattering/boundary value problems (cf.[19, 17, 18]) since we have (1.14) and 
thus z • z — > oo as r — > oo. In the next subsection we study more about Ca{z, 5) by 
carefully checking a proof of (3.7). 

3.2 Asymptotic behaviour 

Define 



A(c,r) = 




The aim of this subsection is to give a proof of the following estimates. 
Proposition 3.1. Let R > and —1<S<0. We have 

\\D^9z * /I|l^(r») < icTXp-'Cs,R\\f\\Ll^,iKr^), |«| < 2, ctX > R. (3.8) 

A change of variables yields 

gz{x) = {cTX{c,T))'^~'^hx{y;uj,uj-^)\y=cTX{c,T)x,x=x{c,T) (3.9) 

where 

K{y- CO, cu^) = -—— / ^ < A < 1. (3.10) 

Given an arbitrary rapidly decreasing function / we give an estimate of hx * f. Using 
a rotation in ^-space, it suffices to consider the case when cu = (0, 1, 0, • • • , 0)^ and a;-*- — 
(1,0, •••,0)^. 

Define 

E = {eeR"]|e + c^^l' = l, 
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E is the set of all real zero points of the complex- valued polynomial 
however, E itself is independent of A. 

Lemma 3.1. Given e > there exists a positive constant Cg independent of A such that 
for all ^ e R'* with c?zsfe(0 ^ ^; /ia?;e 

iMe;A)i>c,ier. 

Proof. Since < A < 1, we have 

IMe;A)|>|Me;i)l- 

Let ICI > 2. We have 

IMe;l)r = (ier + 2a;^-0' + 4(u;-0' 
= + 4(a;^ • Oier + 4(a;^ • 0' + 4(0; ■ 0' 

>ier-4iei^ = ier (i-^) 

Since ^ 7^ and 1)| > for ^ with dists({) > e, the function |Cr^|/i(C; 1)1 is contin- 
uous and positive on the nonempty compact set K^ — {^e R" | |^| < 2 anddist s(0 > 
Thus = inf^gXe 1)1 > ^-^id choosing 

= min ^-^,m^ , 

we obtain the desired estimate. 

□ 

For the treatment of h{^; A) in a neighbourhood of E, we start with the following 
fundamental fact: 

1 1 r e^^'' 



2'Ki{xi + 1x2) (27r)2 J r]i + ir]2 

This yields 



1 1 r e*^-^ 

— d77. (3.11) 



2m{\-^xi + ixs) (27r)2 7 771 + iA-1772 
Lemma 3.2. Let -1 < 5 < and < A < 1. Then 



( f 
Zxf-- ' ^ 



defines a bounded map from L'^_^_^{'R^) to L^(R'^) and its operator norm is bounded from 
above with respect to A. 
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Proof. Prom (3.11) we have 



'Xi + 1X2 . 

where * denotes convolution with respect to variables {xi,X2)- Let g be an arbitrary 
rapidly decreasing function. Since < A < 1 we have 

|A~"^a:i + ix2\ > \x\. 



This gives 



i27r)'\<ZJ,g>\' < 



9{x)f{y) 



A ^{xi - yi) + i{x2 - y2) 



dxdy 



<- (II 



9{x)f{y) 



dxdy 



{xi - yi) + i{x2 - y2) 

Thus everything is reduced to the case when A = 1 and it is nothing but Lemma 3.1 in 
[19]. 

□ 

Having Lemmas 3.1-2 and using a rotation in ^-space, localization and a change of 
variable which are exactly same as Sylvester-Uhlmann's argument [19] for the case z-z = 0, 
we obtain 

WD'^h, * /|U.(R„) < C,||/|U2^^(R„), \a\ < 2. (3.12) 

The constant Cs is independent of A. 

To deduce (3.8) from (3.12) we employ a scahng argument [6]. It follows from (3.9) 
that 

9z*f = {crX)~^ (/iA * /(crA)-i)^^^ , A = A(c, r), (3.13) 

where 

friix) = f{rjx), > 0. 

Let s > 0. We have 

IIA-illLi(Rn) < 7?^+"/'max(l,i?-^)||/|U2(R„), 0<R<r). 
Now from this and (3.13) we have 

\ < 9z* f,9 > \ ^ (ctA)"^| < (hx * /(cTA)-i)crA, 9 > I 
= (ctA)-(2+")| <hx* /(crA)-i,5(crA)-i > | 

< (cr A)"'^^+"^||/lA * /(cTA)-i||L2(R")||fi'(cTA)-i||L2^(R") 

< (CT A)~^^+"^C5||/(c-rA)-i||L2_^^(R«)||5'(crA)-i||L2^(R") 

< (crA)-(^+")C,(crA)^+i+"/^ max(l, i?-(^+i))||/|U.^^(R„)(crA)-^+"/2 max(l, R')\\g\\r^._^^^.) 



{cTX)-'C{S,R)\\f\\L.^^^^nMLlAnr^)- 



This together with the same argument for D°'gz * f yields (3.8). 
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3.3 Uniqueness and Construction of 

Write (1.24) as 

(-A - 2z ■ V)e^ = (ra + fe)e^ - (ra + h) inR". 

Since the solution of (3.4) is unique in L^(R") and has the form (3.6) for / G L^_,_^(R"), 
we have 

= TQz * (oe^) -rg^* a + g^* (be^) -g^^b. (3.14) 

Define 

K,{a) : L^R") 3h^rg,* (ah) e L^(R") 

and 

L,{b) : L^(R") 9 /i ^ * (bh) e L^(R"). 

It follows from (3.7) that both and define bounded hnear operators in L^(R"). 
Rewrite (3.14) as 

(J - Kz{a) - Lz{b))ez = -Tg^ *a- g^^b. 
It follows from (3.8) for \a\ — that 

ll^^(a)^llL2(R") < (cA)"^C5,il||a/i||L2^^(Rn) 

provided cA > i? > 0. Since ||a/i||L2^^(Rn) < || < a; > o||l°°(R")||^||l2(R"), one gets 

||^;j(a)|| < {c\)~^Cs^r\\ <x> a||Lcx,(R„) 

and similarly 

11^^(^)11 < {ctXY^Cs^rW <X> 6||L^(Rn). 

From these we obtain 

||X,(a) + L,(6)|| < (cA)-^C(5,i?)(|| <x>a|Uoo(K.)+T-i <x>6|Uoc(R„)). (3.15) 
Let ci > C2 > and R > 0. Let c and r satisfy 



and 



1 

T > 



^2 



We have c^r > 1 and cA > i?. Now given 0<6<1, R>0 and 77 > choose a large Ci 
in such a way that 



Vyci - ci > C(5,i?)(|| < x> a||Loo(Rn) + clW <x> 6||l-'(R"))- (3.16) 



Since c\> \Jc\ — cl, a combination of (3.15) and (3.16) gives 

\\K,{a) + L,{b)\\<7j. (3.17) 

Choosing a 77 < 1 and 
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for ci satisfying (3.16), specially chosen C2 < ci and R, we obtain the uniqueness and 
existence of the solution of (3.12) and thus those of (1.27), too. 
Sz e LjCR"-) takes the form 

oo 

= - Y,{Kz{a) + Lz{h)f{Tgz *a + gz*b). 

n=0 

It follows from (3.17) and (3.8) for \a\ = that 

oo 

||e^||L2(R„) < ^r/"(cA)~^C(5,i?)(||a||i2^^(R„) + c^||6||i2^^(R„)) 



(3.18) 



n=0 

= Y^{cX)-'C{S, i?)(||a|U.^^(R„) + c^||6|U^^^(Rn)) = 0{{cX)-'). 

Differentiating both sides of (3.14) in the sense of distribution, we have 

(7 - K,{a) - K,{b))dje, = -rg, * dja - K,{dja)e, -g,* djb - L,{djb)e,. 

A similar argument for the derivation of (3.18) yields 

l|9,e,|U2(R„) = 0((cA)-^). 

Continue this procedure for higher order derivatives of in L^(R") and apply the Sobolev 
imbedding theorem to the resulted estimates. Then one concludes that e^ln can be iden- 
tified with a function in C^(Q) and 

||ez||L°°(a) + 1 1 Ve^ 1 1^00(^2) = 0((cA)"^). 

Thus choosing again a large ci, we obtain (1.28). This completes the proof of Theorem 
1.2. 
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